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BASIC STUDIES OF A LOW DENSITY HALL CURRENT ION ACCELERATOR 
by Donald L Chubb and George R. Selkel 

Lewis Research Center ^ 

/ fc ^ 

SUMMARY 

An experimental and theoretical investigation of a low density annular Hall current 
ion accelerator operating with argon gas was carried out. The variation of electron tem- 
perature and number density, ion current density, and plasma potential through the accel- 
erator were experimentally determined as a function of the magnetic field strength. Azi- 
muthal electron Hall current was also measured. Experimentally determined values of 
the axial electron current density parallel to the applied electric field were much larger 
than those predicted by classical electron diffusion, and the ratio of the azimuthal Hall 
current to the axial current was much smaller than expected. As a result, it was con- 
cluded that an anomalous electron diffusion mechanism exists. 

The variation of the electric field with magnetic field could be reasonably predicted 
by an anomalous electron diffusion analysis that made use of the experimentally deter- 
mined value of the Hall parameter (electron cyclotron frequency/electron collision fre- 
quency). The experimental results indicate that energy and momentum transfer to the 
accelerator walls are the primary loss mechanisms for the accelerator. 

INTRODUCTION 

Hall current ion accelerators are currently receiving considerable research atten- 
tion. The use and performance of such accelerators have been discussed by several 
authors for the low density regime (refs. 1 to 4). This regime is defined as that for 
which the mean free path for ion-neutral collisions is much larger than the characteris- 
tic length of the accelerator. Hall accelerators that operate in a higher density range 
where ion-neutral collisions are important, have also been studied (ref. 5). 

The annular Hall ion accelerator to be considered here is illustrated in figure 1. 
Operation of the accelerator with argon gas was investigated. A radial magnetic field 
was produced between an inner iron core and an outer iron ring by a coil wound around 
the inner core. An axial electric field, which accelerates the ions, is established be- 



tween the cathode and anode. The perpendicular electric and magnetic fields establish 
an azimuthal electron current. Since this current is perpendicular to the applied electric 
and magnetic fields, it can be called a Hall current. It is from this that the accelerator 
derives its name. 

Operation of the discharge is maintained by nonequilibrium ionization of the neutral 
argon atoms through collisions with high energy electrons. The electron temperature 
must be high enough to maintain a balance between the ion loss rate and the ion produc- 
tion rate. Ion losses are those to the walls as well as to the accelerated ion beam. The 
elevated temperature is maintained by joule heating of the electron gas. 

Both a theoretical and an experimental study were carried out for this type of annu- 
lar accelerator. The experimental program was directed at understanding the physics 
of the accelerator rather than measuring overall characteristics such as thrust. As a 
result, local macroscopic properties such as number density, electron temperature, and 
plasma potential were measured. From this information conclusions could be drawn 
concerning electron diffusion across the magnetic field and the ion acceleration process. 
Also, the dependence of the accelerator performance on the magnetic field could be es- 
tablished. Probe techniques were employed to determine the plasma properties. From 
the measured properties an approximate acceleration efficiency was computed. 

The theoretical analysis was used to establish the dependence of the electric field on 
the magnetic field. In making the analysis, the effect of inelastic processes (ionization 
and excitation) was considered; however, energy loss to the accelerator walls was not 
considered. 
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THEORY 


Classical Analysis 

In analyzing the annular Hall ion accelerator, a combination of the electron continuity 
and momentum equations was employed. The electron energy equation that was used in- 
cluded ionization and excitation. Such an analysis is similar to the methods commonly 
applied to the positive column in a magnetic field (ref. 6) . 

In appendix B the following equation of motion for electrons is developed (eq. (B15)) 




■ ^^"e + x B) = 


■ ^r.) + “ ^i) 

m + m ®°m ® ° m + m^ ®^m ® 


where the subscripts e, o, and i refer to the electrons, neutrals and ions, respectively, 
and q is the magnitude of the charge of an electron (q^ = -q). (All symbols are defined 
in appendix A. ) The quantities v and v . are the momentum transfer collision 

frequencies for electron- neutral and electron-ion collisions and are given in appendix B 
by equations (BIT) and (B18). 

The following approximations were made in obtaining equation (1): 

(1) The contribution to momentum transfer from inelastic collisions has been ne- 
glected. 

(2) Maxwellian velocity distributions have been assumed for the electrons, ions, and 
neutrals. 



(4) Steady state exists, 3/8t = 0. 

Making the following assumptions 

(1) The accelerator exhibits azimuthal symmetry (3/90 = 0 and = 0) 

(2) Induced magnetic fields are negligible, therefore, only the applied radial field 

exists (B = Bip 

(3) Radial variations in velocity are negligible 

« 

(5) m^ « m^ « m. 

gives equation (1), in scalar form, as 
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Figure 3. - Importance of electron- neutral collisions in argon. 


The results of appendix B were used to calculate v . jv as a function of n /n 

^m/ ® ° 

and T . Data from reference 7 (p. 31) for were used to compute from 

® ^eS. ®°m 

equation (BIT), and the results are presented in figure 2. The results of the calculation 

of V . jv are shown in figure 3. Since In A appearing in equation (B 18) varies 
very slowly with n„ and T^, a constant value of 14 was used, which corresponds to 

c ® 17 ® 

T = 10 K and n. = 10 per cubic meter. For the Hall ion accelerator, the operating 

® ^ 5 o / “2 '“3 

conditions are T^ > 10 K and n^/n^ ~ 10 to 10 . As a result, it can be seen from 

figure 3 that electron-neutral collisions are much more important than electron- ion colli- 
sions. It therefore seems reasonable to neglect terms containing . jv in equa- 

m/ ^m 

tion (2) . If this is done, equations (2b) and (2c) can be written as follows 
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"e y +- 

z \ u. 


^ i s 

^6 rn dz 

® ®°m ®°m ® ® ®°m 


where w = Bq/m is the electron cyclotron frequency. It will now be assumed that 

tJ 


Ue I^QQ ^ « 1 and (1/Wg) ^9Ug Jdzj « 1. Also, it can be shown from experimental 


data that u- v 1 (3u /0z] « 1. As a result, equations (3) and (4) become 


a 1 ^Pe 

u = u 2 E + — ^ u„ s 

®z °z 3z 

e eo eo 6 6 eo 

m m m 


Substituting equation (5) into equation (6) gives the following expressions for u and 


G CO.V, I 0 I 60.V, eel « 

1/2 / m \ 2 

\ eo / \ eo 


H“e»e) \ 
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Tor a successful Hall accelerator, » 1* In this case 


CO V 

e eo 


m 


+ 1 


^e/^eo ^ equations (7) and (8). Also, if (\)/[* "HA 


2 

eo. 


m/j 


and the electron pressure gradient term dn^H^jdz are neglected, equations (7) and (8) 
become: 



e eo 


m 


E_ = 


qB' 



( 9 ) 



( 10 ) 


The azimuthal electron velocity given by equation (10) is just the E x B/B drift velocity 
that the electrons would attain in a collisionless situation. 

Under the same assumptions used to develop the electron equation of motion, the 
following electron energy equation can be obtained 


V 



Vpg + qn^E . Ug 


3m m. 

e-i— (kTj - 

iC m 


(mg+m.)' 


3m m 
e o 

(”e + “o) 


- (kT„ - kT^)v . 

^ m m. 


m / . 


( 11 ) 


where Vj and are the ionization and excitation collision frequencies, and cp^ and 
cp are the ionization and excitation potentials of the neutrals. The development of this 
expression is carried out in appendix C (eq. (C7)) . For the experimental Hall acceler- 
ator, and Tg » Tg,' therefore, the first term on the right side of equa- 

m m 

tion (11) will be neglected, and kT^ will be neglected compared with kT^. Also, it will 

be assumed that the only important energy transfer term is E • u ^ and that 

« E„U- : therefore, equation (11) becomes 
r e.„ z e„ 


®z ^ 



V 

eo eo. 


kT. 


m 


( 12 ) 
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where 


% 


m 


eo 


2m m 
e o 


(m^ + m/ J kT^'-eo * “e) 


(Wl‘'l + El'"ex‘'ex) 


2”e Wl‘'l + E“«’ex''. 


ex 




O 6 60,v, 

Z m 


The quantity k which is a function of T , is the average energy lost by an electron in 
eo V? 

a collision with a neutral. Expressions are developed for Uj and in appendix C and 
are given by equations (C17) and (C18). 

Solving equation (12) for u^ results in 



3 

2 


/c„ kT. 
eo eo^ e 
m 


qE, 


(13) 


The expression that relates electron temperature, magnetic field, and electric field 
is obtained by combining equations (9) and (13) : 



Q ^ 

*5 GO 
2 


kT 


(14) 


The elastic energy loss term 2mg/m^ in the parameter is much smaller than the 
ionization and excitation loss terms for the electron temperatures that exist in the argon 
Hall accelerator. Therefore, equation (14) can be written as 





(15) 


Equation (15) was used to obtain an expression relating E^, B, N^, and T^. By 

choosing a value of T^ and and then determining Uj and (fig. 18, appendix C), 

E can be determined as a function of B. The electron temperature is determined by 
z 

the ionization level necessary to maintain the discharge. It must be maintained at a level 
that gives a balance between the ionization rate and the loss rate of ions and electrons, 
including losses to the walls. Since the simplified analysis presented here does not in- 
clude this wall transfer process, the electron temperature is an unknown parameter in 
determining E„ as a function of B from equation (15). 
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Anomalous Diffusion Analysis 


The electric field strengths predicted by the classical collision theory are much 
larger than those obtained experimentally. Therefore, an anomalous electron diffusion 
across the magnetic field as suggested by Janes and Dotson (ref. 8) and Brockman, Hess, 
and Weinstein (ref. 9) was considered. In the case of anomalous diffusion for 
(^e/^eo following expression for the electron velocity perpendicular to the 

magnetic field, neglecting the electron inertia terms, has been derived in reference 10 
(eq. (33)); 



£L 

B 



9z 


If the pressure gradient term is neglected, equation (16) becomes 



(16) 


(17) 


where a is a constant related to the mean- square deviation of the number density from 
the steady- state value of the number density (ref. 10). Comparison of the classical elec- 
tron velocity equation (eq. (9)) with equation (17) shows that a plays the same role as 
the Hall parameter . As a result, 1/oc might be called an effective Hall pa- 

rameter. 

The electron energy equation is given by equation (12). If the elastic energy transfer 
term 2m /m is neglected in the parameter k , the result for u is 


u 




(18) 


If equations (17) and (18) are combined, the following expression results: 

e2 

+ ( 19 ) 

Comparing equation (19) with its classical counterpart (eq. (15)) shows again that 1/a 
appears in place of the Hall parameter. 


9 



Ceramic center body-^ 


Magnet pole^ 



Figure 4. - Experimental Hall current ion accelerator. 
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Figure 5. - Schematic diagram of experimental accelerator. 






EXPERIMENT 


Apparatus 

A photograph of the experimental accelerator is shown in figure 4. Figure 5 is a 
schematic diagram of the accelerator and associated electrical equipment. The discharge 
was produced in an annular region formed by two concentric tubes. The outer tube was a 
6-inch section of standard 3-inch diameter pyrex tube, and the inner tube was a 3. 8- 
centimeter-diameter ceramic tube that fitted over the inner pole of the magnet. A brass 
flange and gasket provided a vacuum seal at the anode end of the accelerator, and at the 
cathode the 3-inch pyrex tube was connected to a 24-inch-square aluminum box connected 
by a 12-inch valve to the vacuum facility. This vacuum facility is a 15-foot-long 5-foot- 
diameter stainless- steel tank equipped with four 32-inch diffusion pumps. Nominal pump- 

5 

ing speed of the facility is 10 liters per second. 

Argon gas was introduced through four equally spaced holes in the annular stainless- 
steel anode. Gas flow rate was monitored with a rotometer. The pressure of the dis- 
charge was measured at the anode by a McLeod gage. The gage was equipped with a 
liquid- nitrogen- cooled trap. An emitting cathode consisting of twelve 0. 005- inch tanta- 
lum wires was used. The cathode (fig. 4) was constructed much like a bicycle wheel with 
the tantalum wires as the spokes, where the inner negative connection is the hub, and the 
outer positive ring is the rim. All power supplies were standard direct- current supplies, 
except the magnet supply, which had negligible ripple at full 100-ampere output. 

Radial magnetic field strengths in the accelerator were calibrated as a function of 
current in the field coil. A profile of the radial field strength for a current of 9 amperes 

is shown in figure 6. The value of the radial 
field in the middle of the accelerator (200 
gauss) is the value quoted in presenting the 
experimental results. To avoid the hystere- 
sis problem in calibrating the magnetic field, 
the iron magnet was first driven into the 
saturation region and then the calibration 
was made as the current was decreased from 
saturation. Therefore, to use the calibration 
curve, the magnet was first saturated and 
then the current reduced until the desired 
field was obtained. The maximum radial 
field that could be attained in the middle of 
the accelerator was 250 gauss. 

Figure 6. - Radial magnetic field strength profile. 
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Discharge Operation 


The discharge exhibits characteristics similar to those of the positive column of a 
glow discharge (ref. 11, p. 209); that is, an approximately constant electric field exists 
in the discharge, and the plasma is essentially neutral. The discharge was initiated in 

_3 

the pressure range of 10 millimeter of mercury by raising the anode potential to be- 
tween 200 and 300 volts and then switching the magnetic field on and off. Once the dis- 
charge was established, the anode potential and magnetic field were adjusted to the de- 
sired conditions. Before collecting data, the discharge was allowed to reach steady-state 
operation. All the data reported here were obtained at a constant total current of 1 am- 
pere. If it is assumed that the current is uniform across the accelerator, of 1 am- 
pere corresponds to the total current density J,p of 289 amperes per square meter. All 
potentials reported were measured with respect to ground. 
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(c) Button probe and circuit 
Figure 7. - Diagnostic probes. 






• Operation of the discharge is strongly dependent on the cathode temperature, indi- 
cating that the emitted cathode current is emission limited rather than space- charge 
limited. By raising the cathode temperature, the total current increased while the anode 
potential decreased. In obtaining the data presented here, the cathode temperature was 
maintained at a constant value throughout the interval when data were being recorded. 


DIAGNOSTIC METHODS AND EQUIPMENT 
Langmuir and Emitting Probes 

Measurements of electron temperature, electron number density, and plasma poten- 
tial were made with a 0. 001-inch-diameter cylindrical tungsten Langmuir probe and a 
0. 00 5- inch- diameter tungsten emitting probe. Both probes are illustrated in figure 7. 
Operation of an emitting probe is described in reference 12. As the probe is heated, the 
floating potential rises to a plateau. At this point, the emission current from the probe 
equals the random electron current arriving at the probe. As a result, the floating poten- 
tial of the probe closely approximates the plasma potential. The emitting probe was 
used only in obtaining the plasma potential. 

It was assumed that the magnetic field did not affect the Langmuir probe operation. 
This appears to be a reasonable assumption since the measured Hall current through the 

radial cross section A^ of the volume, 
from which the probe could easily draw 
current, was large compared with the 
measured saturation electron current 
(fig. 8). 

The Langmuir probe data were re- 
duced by two different methods: first, 
the relation between the electron cur- 
rent and the electron distribution func- 
tion derived by Druyvesteyn (ref. 13), 
and second, a method described subse- 
quently as "corrected Maxwellian. " 

In determining the electron number 
density n^ by either method, the 
plasma potential used in the calculation 
was obtained with the emitting probe. 
Associated with the anomalous 
electron diffusion across the magnetic 


r 
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field are large oscillations in the 
plasma potential. For low magnetic 
fields (B < 50 gauss), the oscilla- 
tions are small. As pointed out in 
reference 14, oscillations near the 
plasma potential will cause a reduc- 
tion in the electron saturation cur- 
rent measured by a Langmuir probe. 
In figure 9, typical Langmuir probe 
curves are plotted for magnetic 
fields of 50 and 250 gauss. The 
plasma potentials VpL shown in 
figure 9 were determined by the 
emitting probe. If the electron dis- 
tribution function is Maxwellian, the 
probe curve should increase linearly 
to VpL and then level off (ref. 15). 
As shown in figure 9, however, the 
probe curve has a break at a lower 
value Vp, rather than VpL- This 
result can be interpreted in two 
ways. The most obvious conclusion 

is that the distribution function is not Maxwellian. Secondly, it could be concluded that 
the oscillations in the plasma potential cause a premature break in the probe curve, but 
that the distribution function is still Maxwellian. In this case, the true electron satura- 
tion current can be taken as i^ (fig. 9). If the probe curve is interpreted in this manner, 
the electron number density and temperature can be determined as usual (ref. 15). The 
probe data reduced this way is referred to in the figures as corrected Maxwellian. 

If the effect of the plasma oscillations is ignored, the data can be reduced by the 
Druyvesteyn method in the following manner. If an isotropic distribution function is as- 
sumed to exist for the electrons, the second derivative of the curve for electron current 
against probe voltage will be proportional to the electron distribution function (ref. 13) 



Figure 9. - Typical Langmuir probe data. 


n(p(U)=-^ 
® Aq 


/2m U 
e 


. 1/2 


d 1. 


( 20 ) 


dv: 


where U = V^it - Vt^ is the difference between the plasma potential and the probe po- 

rL, P 2/2 

tential, A is the probe surface area, d ig/^Vp is the second derivative of the electron 
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* current collected by the probe, and ng(/>(U) is the electron density within the energy in- 
■ crement dU. 

By definition, the electron number density is 


“e' jf” 


Substituting equation (20) into equation (21) and integrating by parts yield 


( 21 ) 



and the first term in the parentheses vanishes because (dig/dU)-jj_^ = 0. Therefore, the 
electron number density is 



(23) 


The electron current i^ was taken as (i - i.), where i is the measured probe current 
and i. was estimated graphically by drawing a straight line through the point determined 
by zero current and plasma potential and tangent to the ion saturation portion of the probe 
trace. To reduce this equation (23) further, a relation must be obtained between i^ 
and U. For all experimental probe traces that were recorded in cases where B > 0, 
this relation could be well represented by two exponential curves. (Fig. 9 is a typical 
plot of the logarithm i^ against the probe potential. ) Therefore, the electron current to 
the probe can be approximated as follows 


-vU 

= V 

0 < U < U* 

= ile-SU 

s 

U*< U < « 


(24) 


where U* = Vpj__ - Vp, 6 and y are the slopes of the two straight-line approximations 
to log ig, ig is taken as the electron current at the plasma potential VpL’ and ig is 
the electron saturation current obtained by extrapolating the first part of the probe trace 
to the plasma potential VpL- When B « 0, the plasma potential determined by the break 
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in the probe trace is slightly higher than the plasma potential measured with the emitting 
probe. Using the plasma potential determined by the emitting probe for cases where 
B 0 would therefore appear to be a conservative estimate. 

Using equation (24) in the number density expression (eq. (23)) results in 



where the integrals are incomplete gamma functions (tabulated in ref. 16, p. 978). 

The mean random electron energy, or electron temperature, can also be determined 
from the probe trace, since by definition 



(26) 


where (kT /q) is the experimentally determined electron temperature. This temper- 
ature is not necessarily equivalent to the electron temperature used in the theory, since 
in deriving the expressions for v , Uj, and it was assumed that the distribution 
function was Maxwellian. 

Substituting equation (20) into (26) and integrating by parts give 


3 

2 




where the condition (dig/dU)y_^ = 0 has been used. Finally, use of equation (24) in 
equation (27) yields 






l/2e-x 


dx 




xl/2e-x 



(28) 


where the integrals are the same as in equation (25). 
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Hall Current Measurement 


As a check on the methods described for measuring the electron number density, a 
method based on knowing the electron Hall current in the accelerator was also used. The 
azimuthal Hall current was measured in a way similar to that used in reference 9. When 
the discharge is turned off, the collapsing azimuthal Hall current induces a voltage in a 
coil wrapped around the outer pyrex tube. To pick up this induced voltage, a 1000-turn 
search coil was wrapped around the pyrex tube just downstream of the outer pole of the 
magnet. The coil was located at an axial station, z = 5. 08 centimeters (fig. 5). The sig- 
nal from the coil was fed into an integrating circuit, and the result was displayed on an 
oscilloscope. The search coil was calibrated by simulating the Hall current with a coil 
placed between the poles of the accelerator magnet. By switching off a known current in 
the calibration coil, the recorded oscilloscope reading was plotted as a function of the 
simulated Hall current. The steady magnetic field produced by the iron magnet of the ac- 
celerator had a negligible effect on the calibration. The effect of a nonuniform Hall cur- 
rent was simulated by moving the calibration coil. By placing this coil at several axial 
positions between the magnet poles with the same value of current and comparing the os- 
cilloscope readings, the effect of nonuniformities in the azimuthal Hall current could be 
determined. A maximum change of 20 percent from the reading when the calibration coil 
was centered under the magnet poles resulted when the coil was moved to a position with 
one edge at the cathode position. This would simulate the extreme case of having most of 
the Hall current occurring between the cathode and the downstream edge of the magnet. 

The measured electron Hall currents were used to obtain a mean number density in 
the following manner. From equation (10), which neglects the electron pressure gradient 
and assumes ^ 





B 


Therefore, the mean number density can be obtained from measured mean values of 
^Hall’ ®z’ 


B 

To convert the measured Hall current into a current density, it was divided by the area 
bounded by the inner and outer tubes, the anode, and the downstream edge of the magnet. 
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This area equals 9. 67 square centimeters. 

If the electron Hall current density and the axial electron current density are known, 
the mean effective Hall parameter can be computed. The following expression is 

obtained from equation (5) : 



= 

e e 


Button Probe 

Axial electron current densities were determined by measuring the axial ion current 
density and then subtracting this result from the total current density. The ion current 
density was measured with a 0. 5-centimeter-diameter guard- ringed tungsten button probe 
shown in figure 7(c) (p. 12). Biasing the button negative with respect to the plasma poten- 
tial results in ions being collected by the button, while electrons are repelled. To com- 
plete the electrical circuit, an electron collector is placed near the button. As the button 
becomes more and more negative, eventually all the electrons will be repelled, and only 
ions will be collected. There will be some value of the probe potential beyond which no 
further ion current will be collected by the button. The value of the current at this point 
is an upper bound on the axial ion current moving through the plasma. In operating the 
probe, both the button and the button guard ring were maintained at the same potential. 
When this was done, the probe current saturated, indicating that any additional ion cur- 
rent attracted by the sheath growth was collected by the guard ring. 


Neutral Number Density Measurement 

The McLeod gage pressure measurements were used to establish the neutral number 
density. The neutral atoms are essentially in free molecular flow throughout the dis- 
charge; examination of similar cold flows indicates a small (30 percent) pressure drop 
through the accelerator. This coupled with the knowledge that the discharge is only 
slightly ionized indicates that the neutral density is approximately constant. Since the 
mean free paths of the particles in the accelerator are larger than the characteristic 
length of the apparatus, the condition of equal pressures in the accelerator and McLeod 
gage cannot be used in determining the neutral density; however, conservation of mass 
must exist. If it is assumed that the electrons and ions have completely recombined in 
passing through the tube connecting the accelerator and the McLeod gage, the conserva- 
tion of mass can be expressed as follows 
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( 29 ) 


* P*v’=DV +P.V. + PV 

o e ^11 o 

where the primed quantities indicate conditions in the McLeod gage and the unprimed 
quantities indicate conditions in the accelerator and = n^m^, = n^m^, p^ = N^m^. 

If equal ion and electron densities and m^« m^ are assumed, the neutral number density 
is obtained from equation (29) 


No 



(30) 


v' 


For free molecular flow, v = VskT/rm. This free molecular flow speed was used for 
. ^ and Vq in equation (30) . For the ion velocity required for a stable sheath, given 
by the expression v. = ^kT^/m^ (ref. 17), was used. Since current conservation 
requires equal arrival rates for electrons and ions (riv^ = riv^) and m^ is neglected com- 
pared with m^, equation (30) becomes 


No 



( 31 ) 


Typically, ~ 10 per cubic meter and ny^Tg/ST^ ~ 10 per cubic meter; there- 
fore, the second term in the numerator of equation (32) can be neglected, and the neutral 
density in terms of the measured McLeod pressure becomes 


= N’ 
o o 



; Vt’T^ 

f o o 


(32) 


Before can be computed, the temperatures T^ and T^ must be known. The tem- 
perature in the McLeod gage should be nearly room temperature; therefore, a value of 

300° K was used for T* . The neutral temperature in the discharge T is not known, 

O Q ^ o 

but in order to compute N^, a value of 500 K was chosen for T^. This value of 500 K 

should be higher than what actually exists, therefore, the computed neutral number den- 
sities might be slightly low. 
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RESULTS AND DISCUSSION 


Electron Distribution Function 

If a Maxwellian distribution exists for the electrons and no plasma oscillations are 

present, the Langmuir probe curve, log i^ against Vp, would lie on a straight line for 

Vp<VpL- As shown in figure 9 (p. 14), however, log i^ is not a single linear function, 

but can be closely approximated by two linear functions. It might therefore be concluded 

either that the electron distribution function is not Maxwellian for B > 0 or that plasma 

oscillations exist. More likely, a combination of these two conditions exists. In other 

words, the distribution function is not Maxwellian, and the probe curve exhibits a reduced 

electron saturation current (i in fig. 9) due to plasma oscillations. 

s 

Since a Maxwellian distribution was assumed in deriving the expressions for v , 

^°m 

Vj, and V , the theoretical analysis of the accelerator will be in error. The high-energy 
tail of the distribution is the region where the largest departures from a Maxwellian are 
likely to exist. As a result, the ionization and excitation collision frequencies, which de- 
pend primarily on the high-energy tail of the distribution, will show the greatest error be- 
cause a Maxwellian distribution was assumed. If the high-energy tail of the distribution 
function is below that of a Maxwellian distribution, the values of and Uf based on a 

0X 1 



0 .2 .4 .6 .8 1.0 1.2 0 .2 .4 .6 .8 1.0 1.2 


Dimensionless axial distance, z/L 

(a) Plasma potential. (b) Electron temperature. 

Figure 10. - Variation of plasma properties through accelerator. Total applied voltage, 162 volts; total current density, 289 amperes per 
square meter; magnetic field strength, 200 gauss; neutral number density, 5. 2xl0l9 per cubic meter; mass flow rate, 10‘^ gram per 
second; accelerator length, 7. 62 centimeters. 
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Maxwellian distribution will be larger than the values obtained for and u, with the 
' correct distribution function. 


Experimental Variation of Plasma Properties Through Accelerator 

In figure 10 the experimentally determined variations of plasma potential, electron 
temperature, number density, ion current density, and axial ion velocity through the ac- 
celerator are presented. The electron number density and temperature were determined 
by both the Druyvesteyn and corrected Maxwellian methods. It can be seen from figure 
10(a) that the plasma potential varies almost linearly with axial position through much of 

the accelerator and as a result, E_ is approximately a constant. A 50- to 60-volt sheath 

z 

was observed to exist at the cathode for all magnetic fields investigated. 

As shown in figure 10(b) the electron temperature is a generally increasing function 
in moving from the cathode to the anode. The corrected Maxwellian electron tempera- 
tures are lower than those determined by the Druyvesteyn method. For z/L > 0. 7 the 
magnetic field is small since these stations are outside the magnet poles, and the 
Druyvesteyn data and corrected Maxwellian data tend to approach agreement. 

Figure 10(c) shows that the electron number density increases rapidly through the 
first one-third of the accelerator and then decreases up to the cathode. The initial in- 



Dimensionless axial distance, z/L 

(c) Electron number density. (d) Ion current density. Average ion current density, 

82. 8 amperes per square meter. 

Figure 10. - Continued. 
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crease in number density, despite the fact that ions are being accelerated, indicates that 
ionization is predominant for the first one-third of the accelerator. Near the peak, the 
corrected Maxwellian electron densities are nearly an order of magnitude greater than 
Druyvesteyn values. Just as in the case of temperature, the Druyvesteyn and corrected 
Maxwellian densities approach closer agreement beyond the magnet poles (z/L > 0.7). 

The results of figures 10(b) and (c) show the importance of the electron pressure gra- 
dient term, which was neglected in the theoretical analysis. For z/L < 0. 2, both the 
Druyvesteyn and the corrected Maxwellian data show that 




800 V/m 


compared with E = 1640 volts per meter. However, for z/L > 0. 2, both sets of data 
z 

show that (1/qn^) [^(ngkTgySz] < 400 volts per meter. Neglecting this term compared 
with E therefore seems justifiable for z/L > 0. 2, but not for z/L <0.2. 

Z 

The ion current density measured with the button probe is shown in figure 10(d). 
Again, the rising ion current density indicates the importance of ionization for the first 
one-third of the accelerator. For 0. 5 < z/L < 0. 8, the ion current density decreases. 


The decrease observed in the region 0. 5 < z 
to the walls rather than recombination. The 



0 .2 .4 .6 .8 1.0 1.2 


Dimensionless axial distance, z/L 
(el Axial ion velocity. 

Figure 10. - Concluded. 


/L < 0. 8 is probably caused by radial flow 
recombination rate is negligible for the 
electron temperatures and number densities 
existing in the argon Hall accelerator 
(ref. 18). 

With the assumption of equal charge 
densities, the results of figure 10(d), and 
the number density from figure 10(c), the 
ion velocity was computed and is shown in 
figure 10(e). The criterion for a stable 
sheath (ref. 17) requires a directed ion ve- 
locity toward the button probe of approxi- 
mately y kTg/m. . Axial ion velocities de- 
termined from this criterion for both the 
corrected Maxwellian and Druyvesteyn data 
are between 4 to 7x10 meters per second. 
Such velocities would result merely by the 
presence of the probe even if the plasma 
were quiescent. This does not mean that 
the probe caused the measured ion veloci- 
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ties, but it does determine a minimum ion velocity that can be measured. It is evident 
from figure 10(e) that the ions are essentially unaccelerated. There are two likely ex- 
planations for such a result. First, charge exchange may be responsible. The charge- 
exchange cross section data for argon (ref. 19, p. 164) were used to compute a mean 

free path for charge exchange. For all ion energies greater than 25 electron volts, the 

- 1 fi 

argon charge- exchange cross section is approximately 30x10" square centimeter. 

19 

With this value of the cross section and a neutral density of 5. 2x10 per cubic meter, 
the following mean free path results: 



(5. 2) (30) 


It can be concluded from a comparison of this value of X with the accelerator length, 

L = 7. 62 centimeters, that charge exchange is an important mechanism for limiting the 
ion velocity. Charge exchange is not necessarily a serious inefficiency, since the high- 
speed neutrals produced are just as efficient in producing thrust as ions of equal speed. 
The inefficiency that does result from charge exchange is the nonuniformity in velocities 
at the accelerator exit. 

The second possible explanation for the ion velocity results shown in figure 10(e) is 
that the ions are being accelerated to the accelerator walls where they recombine. Since 
the radial magnetic field has no effect on the radial motion of the electrons, their centrif- 
ugal acceleration will force them to the outer wall of the accelerator. In addition, a ra- 
dial ambipolar electric field will be established that will also accelerate ions to the walls. 
With the assumption of a Boltzmann distribution for the electron density, Langmuir and 
Tonks (ref. 20) have shown, by solving Poisson’s equation including ion production, that 
a potential drop of about kT^/q to 2kTg/q can be established in a cylindrical plasma. 
This potential drop will accelerate ions to the cylinder walls. 

Janes and Dotson (ref. 8) have obtained experimental results for the variation of the 
ion current and electron number density through a similar accelerator. They measured 
the ion current J. with a plane double Lai^muir probe and interpreted the measured 

z 2 

force acting on a small quartz plate as the momentum of the accelerated ions m^n-u^^ . 

z 

When equal charge densities were assumed and the ion current measurement and the ion 
momentum measurement were combined, the electron number density n^ and the ion ve- 
locity u. were obtained. The results of reference 8 are reproduced in figure 11. If 
^z 

charge exchange and location of ionization are important, the interpretation of results 

will be affected since the momentum measurement will also include the neutral momen- 
2 

turn m N u and the force resulting from the ion velocity distribution, which were as- 
o o o_ ’ 

Z 

sumed to be negligible. The results in figure 11 for n^ and J. variations with axial 
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Figure 11. - Experimental data of reference 8. Mass flow rate, 0. 002 gram per second; total current, 
10 amperes; total applied voltage, 400 volts. 


distance are very similar in shape to those of figures 10(c) and (d); however, the results 

for u. are quite different from those shown in figure 10(e). They indicate that u. is 
h z 

a monotonically increasing function from the anode to the cathode. The value of u. at 

the cathode is given as y^qV^p/m., which would be true only if the ions were accelerated 

through the potential V,p. Obviously, if the location of ionization and charge exchange 

were considered the results presented for u^ would be reduced. If charge exchange 

z 

were to predominate, it would be expected that u^ should decrease near the cathode. 

z 

An interesting conclusion can be made from the force measurement of reference 8. 

2 

These results, shown in figure 11 as m.n.u. , imply that the force increases at the be- 

ginning half of the accelerator and then decreases through the remainder of the acceler- 
ator. Since an ion accelerating electric field is maintained through the entire acceler- 
ator, it seems reasonable to expect the force to be an increasing function, even if charge 
exchange is important, unless the plasma momentum is going to the walls. Because 
Janes and Dotson’s results do not show an increasing force, it can be concluded that the 
plasma momentum is being lost to the walls. This conclusion was also shown from the 
results shown in figure 10(e) . 

The importance of the electron acceleration term Su^ hz can be evaluated by using 
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dafa presented in figures 10(c) and (e). From the definition of the total current, an ex- 
pression for the electron velocity can be derived if equal charge densities are assumed: 


Therefore, 


and 






j 

_Jz ^ 

dz q 0z 


(33) 


(34) 


(35) 


From the results of figures 10(c) and (e), the first term in equation (35) can be shown to 

be small compared with the second term. Therefore, the number density curves (fig. 

10(c)) show that the largest values of 0u /dz occur in the region z/L < 0. 2. In this 

z/ 

region 

0U_ T 

J. is- 10« 

0z q ^2 0z 

In order to neglect the term 0u^ jdz in equation (4), the following condition must be 
satisfied: ^ 


0U, 


V.. dz 
m 


« 1 


For 


IQ 

Nq ~ 5x10 per cubic meter and kT^/q > 10 electron volts figure 2 shows that 

1 ^ I -f /-V-M y-v 


eo 

m 


10 . Therefore, 


du 


~ = 0 . 1 


^eo^ 10 ^ 
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As a result, it seems reasonable to neglect the term 



in equation (4). 


Experimental Evidence of Anomalous Diffusion 

To check the validity of the classical expression for electron current density across 
the magnetic field (eq. (7)), the data of figure 10 were used. If (^e/^eoj » 1 

neutral velocity term is neglected, equation (7) gives the following result for the electron 
current density: 



"^e^e'^eo 


•^Ve = 


m 


B 








dz 


(36) 


Both Druyvesteyn and corrected Maxwellian data were used to compute J . The largest 

values of J occur at z/L = 0. 3 for the Druyvesteyn data and at z/L = 0. 35 for the 
®z 

corrected Maxwellian data. At these axial stations, the electron current density from 
Druyvesteyn data at a magnetic field of 200 gauss is 

/"j ^ =2.23A/m^ 

and from corrected Maxwellian data is 

/jg \ = 21. 2 A/m^ 

^ ^^MAX 


In computing J , the values of were obtained from figure 2 (p. 4) . It should be 

z m 

remembered that v is based on a Maxwellian distribution. However, unless the 
®°m 

distribution function is far from being Maxwellian, the correct value of the collision fre- 
quency, which is obtained by integrating over all velocity space, should not be much dif- 
ferent from one based on a Maxwellian distribution. 

To obtain the experimental value of J it was assumed that the total current was 

®z 

uniform across the accelerator, therefore, Jrp = 289 amperes per square meter. Thus, 
from figure 10(d) at z/L = 0. 35, the experimentally determined electron current densi- 
ties are as follows: 
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at z/L = 0. 3 


(jg \ = 179 A/m2 

^'^exp 

^Jg j = 175 A/m^ at z/L = 0. 35 
^ exp 


Based on the large differences between the calculated current densities and the experi- 
mental values, it would appear that there is an anomalous mechanism which allows the 
electrons to diffuse across the magnetic field at a rate much greater than that predicted 
classically. 

Further evidence leading to the conclusion that an anomalous electron diffusion ex- 
ists was obtained by measuring the effective Hall parameter The ion current den- 

sity distribution was measured for several values of the magnetic field with the button 
probe. In all cases the distribution had a form similar to that shown in figure 10(d). 

From these distributions, average ion current densities J. were computed and the re- 

^z 

suits combined with the measured Hall currents to compute an average value of ^gTg 
for the accelerator 


= 

e e 


^6 


'Hall 


•^e„ '^T ■ '^i 

z z 


These values of co are listed in table I. 
e e 

An effective value of can also be obtained from the measured plasma proper- 

ties. Neglecting the pressure gradient term in equation (36) gives the following result: 


qn^ E„ qn^ 

, , ^ ^ e z ^ e z 


(37) 


Experimentally determined values of n and J. were used in equation (37) to compute 

® ^z 

the effective Hall parameters at z/L = 0. 367 for several magnetic field strengths. The 

results from both Maxwellian and Druyvesteyn reduced data are shown in table I. For 

19 

all cases listed, J.p = 289 amperes per square meter and = 5. 2x10 per cubic 
meter . 

The values of a for the Druyvesteyn data and 'a from the Hall current data com- 
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Figure 12. - Dependence of electron number density on magnetic field. Total cur- 
rent density, 289 amperes per square meter; neutral number density, 5. 2x10^^ 
per cubic meter; accelerator length, 7. 62 centimeters. 


pare closely with the value of a = 1/3 presented in reference 8. If figure 2 is used to 
obtain u , the classical value of ^af^eo should be of the order of 100 for the elec- 

tron temperatures, neutral densities, and magnetic fields obtained in the argon Hall ac- 
celerator. As a result of the large discrepancies between the measured values of 

and J and classically computed values of these parameters, it appears that some 
z 

anomalous phenomena occur to produce enhanced diffusion of the electrons across the 
magnetic field. 


Dependence of Plasma Properties on Magnetic Field 

Experimental . - In figure 12 the experimentally determined variation of the electron 
number density with magnetic field for constant total current density J,p = 289 amperes 
per square meter at axial stations z/L = 0. 033, 0. 367, and 1. 03 is presented. Only the 
Druyvesteyn determined values for the density are shown. Maxwellian values for the 
density are nearly an order of magnitude greater, but show the same variation with mag- 
netic field. Also shown in figure 12 is the mean number density determined from the 
measured Hall current. As can be seen, there is good agreement between the number 
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Magnetic field strength, B, gauss 

Figure 13. - Experimental electron temperatures. Total current density, 289 am- 
peres per square meter; neutral number density, 5. 2x10^^ per cubic meter. 

densities determined by the Druyvesteyn method and the Hall current measurement. 

The experimentally determined dependence of the electron temperature and the elec- 
tric field on magnetic field is illustrated in figures 13 and 14, respectively. Both 
Maxwellian and Druyvesteyn data show some increase in electron temperature with mag- 
netic field. From figure 14 it appears that the electric field increases with the magnetic 
field. 

Theoretical . - Results for the electric field that were obtained from the classical 
diffusion analysis (eq. (15)), are shown in figure 15. Comparison of the experimental 
results in figure 14 with those in figure 15 shows that the results for classical diffusion 
predict much larger electric fields than those measured. 

Theoretical results for the electric field under the assumption of anomalous electron 
diffusion are shown in figure 14. The electric field was calculated from equation (19) 
with a = 1/3 and = 5. 2x10^^ per cubic meter. Figure 14 shows that an electron tem- 
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Magnetic field strength, B, gauss 

Figure K - Comparison of experimental and theoretical electric field for anoma- 
lous diffusion. Effective Hall parameter, 3,- neutral number density, 5. 2xl0W 
per cubic meter. 


Electron tern 
perature. 


Q do 60 80 100 

Magnetic field strength, B, gauss 

Figure 15. - Variation of electric field with mag- 
netic field for classical diffusion. 


perature between 10 and 15 electron volts gives 
good agreement between the experimental elec- 
tric fields and the fields calculated. From fig- 
ure 13 it is evident that these electron tempera- 
tures are generally in agreement although slight 
ly less than Druyvesteyn temperatures and do 
not agree with the corrected Maxwellian temper- 
atures. Part of the discrepancy for the Druyves 
tian case may result from the derivation of the 
ionization and excitation collision frequencies 
and As pointed out previously, and 

^ex primarily on the high-energy tail of 

the distribution function. If the actual distribu- 
tion is depleted in this region, the actual tem- 
perature must be higher than a Maxwellian tem- 
perature to produce the same value of and 
fgjj. Therefore, the temperature will also be 
higher to produce the same electric field as a 
Maxwellian distribution (eq. (19)). 







40 80 120 160 200 240 280 

Magnetic field strength, B, gauss 

Figure 16. - Experimental accelerator efficiency. Total current den- 
sity, 289 amperes per square meter; neutral number density, 5. 2x10^^ 
per cubic meter. 


Experimental Accelerator Efficiency 


From the experimentally determined values of the ion current density and the 
Druyvesteyn determined electron number density at z/L =1.0, an accelerator efficiency 
can be computed. To make the calculation, it was assumed that, at the accelerator exit, 
the plasma has equal charge densities n^ = n^ and that the ion velocity and number den- 
sity were uniform. It was also assumed that the total current density Jrp was uniform 
across the accelerator. Under these conditions, the efficiency is 


V 


exp 


^ r ivL 


/z=L 




rv. 

Ir, T Irr T 

1 z=L z=L 
2q 


(38) 


In figure 16, the experimental efficiency is presented. It can be seen that the magnetic 
field had little effect on the efficiency for B > 100 gauss. This seems to substantiate 
that momentum and energy transfer to the walls are the important loss mechanisms. 

With the radial magnetic field, both electrons and ions are free to move along the field 
lines to the walls, so that the magnetic field strength should have little effect on energy 
and momentum transfer to the walls. It should be remembered that the efficiency defined 
by equation (38) does not include the momentum of high-speed neutrals produced by charge 
exchange. If this were included the efficiency would be higher. 


CONCLUSION 

The experimentally determined axial electron current moving perpendicular to the 
magnetic field lines was at least an order of magnitude greater than the classically com- 
puted result. This result, together with the fact that measured values of the effective 
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Hajl parameter were much lower than expected, implies that an anomalous electron diffu- 
sion across the magnetic field exists. Theoretical results for the variation of electric 
field with magnetic field that include anomalous electron diffusion agree well with the ex- 
perimental results. At constant current it was found that the electric field increases 
with magnetic field. 

Experimental results indicate that the ions are not being accelerated to velocities 
that would be attained if the ions were to undergo electrostatic acceleration through the 
applied potential. Charge exchange is a possible explanation for this result. However, 
since Janes and Dotson measured a decrease in the axial plasma momentum, it appears 
that the plasma momentum is being lost to the walls. 

The experimentally determined accelerator efficiency is essentially independent of 
magnetic field equal to or greater than 100 gauss. The independence substantiates the 
conclusion that ion momentum and energy transfer to the walls are the most important 
loss mechanisms. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 17, 1965. 
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SYMBOLS 


probe surface area 

magnetic field strength 

particle random velocity, v - u 

electric field strength 

velocity distribution function 

relative velocity between two 
colliding particles 

total current 

electron current collected by 
Langmuir probe 

current collected by Langmuir 
probe 

current density 

total current density 

Boltzmann constant, 

1. 3805x10"^^ J/°K 

accelerator length, 7.62 cm 

particle mass 

neutral number density 

number density 

1 /* 2 3 

pressure, — m / c F d v 

3 *'v 

heat flow vector, 

— mf c^ c^ F d^v 
2 *'v 

electron charge, 

1.602x10“^^ C 

inner radius of accelerator 


T 

U 


u 

V 

V 


V 

z 

a 

A 

€ 


V 

e 



outer radius of accelerator 
radial coordinate 
Debye length 
temperature 

difference between plasma poten- 
tial and Langmuir probe poten- 
tial, VpL - Vp 

1 /* — 3 

average velocity, - / v F d v 
n *'v 

voltage 

center -of- mass velocity in binary 
collision 

total applied voltage measured 
with respect to ground 

total particle velocity, u + c^ 

axial coordinate 

anamolous diffusion parameter 

inelastic energy loss increment 

azimuthal scattering angle 

accelerator efficiency 

azimuthal coordinate 

average energy lost by electron 
in colliding with neutral, 
eq. (12) 

reduced mass 

electron- ion momentum transfer 
collision frequency 





t'eo electron-neutral momentum 

transfer collision frequency 

t'ex excitation collision frequency 

i/j ionization collision frequency 

cr^ differential cross section 

CTm total momentum transfer cross 

section 

CTrj, total cross section 

T time 

T„ effective electron collision time 

e 

^ex excitation potential for first ex- 

cited level of neutral 

cpj ionization potential 

X scattering angle 

stress tensor, m fccF d^v 
Jy 

cOg electron cyclotron frequency 

Subscripts: 

coll collision 


DR 

Dru 3 Tvesteyn data 

e 

electrons 

e£ 

elastic process 

ex 

excitation process 

exp 

experimental 

Hall 

Hall current 

I 

ionization process 

i 

ions 

MAX 

corrected Maxwellian data 

0 

neutrals 

P 

probe 

PL 

plasma 

r 

radial direction 

s 

species s 

t 

species t 

z 

axial direction 

e 

azimuthal direction 
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APPENDIX B 


ELECTRON EQUATION OF MOTION^ 

K the zeroth moment of the Boltzmann equation is taken the following result is ob- 
tained: 


at 








(Bl) 


where the right side of equation (Bl) represents the rate of production of species s by 
collisional processes such as ionization, and F_(v„, r, t) is the distribution function of 

o o 

species s. If no production occurred, the right side would be zero, and the usual conti- 
nuity equation would result (eq. (6-8) of ref. 7, p. 157). Other terms appearing in equa- 
tion (Bl) are the number density n^, the particle mass m^, the particle velocity v^, 
and the average velocity u„. 

O 

Taking the first moment of the Boltzmann equation, as in reference 21 (eq. (6-9), p. 
157), results in the following expression: 


0(n m u ) ^ 


at 


+ V . iTg - ngqg(E + u g X B) 



(B2) 


The term on the right side of equation (B2) represents the momentum gained by species s 
in collisions with all the other species t in the plasma, and ^s^s^s ^ 

'^s 

the stress tensor of species s. Also appearing in equation (B2) are the charge of species 
s, qg, electric field E, and magnetic field B. 

If equation (Bl) is multiplied by Ug and the result subtracted from equation (B2), the 

^The development of the collision integrals in appendixes A and B follows the work 
of Dr. Frederic A. Lyman presented in unpublished lecture notes at Lewis. In his work, 
however, only elastic collisions are considered. 
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following equation of motion for species s is obtained 


0U 


s 


iUgns — ^ + nging(Ug • V)Ug + V ■ - iigqg(E + x B) = 

9t 


■u 


0T 


- 

d 

s s s 


‘'coll 


(B3) 


where the definition v „ = u „ + c_ has been used in equation (B2). In order to obtain a 

s s s 

useful form of equation (B3), the right side of equation (B3), which represents the mo- 
mentum gained by species s in collisions with the other species in the plasma, must be 
evaluated. 

From reference 22 (p. 62, section 3. 51), the rate of change in momentum of species 
s in collisions with all other species can be written as 




(B4) 


where a prime denotes the condition after a collision, g is the relative speed |v^, - Vg|, 

X is the scattering angle, e is the azimuthal scattering angle, and is the differential 
cross section. The following notation is used for integration over velocity space: 


/.+00 ^+oo ^+oo 

*■ • L L 

The first step in carrying out the integration of the right side of equation (B4) is re- 
lating (ic ' - c ) to g, X) For a binary inelastic collision, where any particles 

s s 

that may be produced are assumed to have negligible momentum and energy, the conser- 
vation of momentum and energy can be expressed as follows: 









+ 1 ^ 


where A is the eiwrgy lo^t in the inelastic process. If the center-of-mass velocity is 

- ™s^s + *”t^t - 

defined as V = — , the conservation of momentum shows that V = V . Writing 

-f ^ _fj*s + "^t 

^s’ ^s’ ^t’ ^t terms of the center-of-mass velocity and the relative velocities 

results in 


• -"^t- 

v_ = V - -i g 

® M 


V = V g 

® M 


Vx = V+— 5g 
^ M 


V ’ = V + -^ g* 
^ M 


(B5a) 


(B5b) 


(B5c) 


(B5d) 


where M = m^ + m^ and g = v^ - v^. Substituting equations (B5) into the energy equation 
results in the expression 


= g 



(B6) 


where /i = (mgmj.)/(mg + m^ is the reduced mass. Also, from equations (B5a) and (B5b) 
the following relation is obtained: 


m. 

V* - V =— ^ (g - g’) 

® ® M 


Since v„ = u _ + c„ and v * = u„ + the following expression is obtained; 

s s s s s s 



m. 

C = — (g - g') 
® M 
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where £ and m are unit vectors. After substitution of equation (B7) into equation (B4), 
integration over e , if a is assumed to be independent of e , eliminates the I and in 
components. As a result, 



where the integral on the right side is the total momentum transfer cross section (t^^. 

For elastic collisions, A = 0, and the elastic momentum transfer cross section is 

^7T 

(s) = 2 t 7 / (l-cosx)cT. (g,x)sinxdx (B8) 

e£ Jq %IL 

where Uj is the differential scattering cross section for elastic collisions. 

%St 

For the inelastic processes, ionization and excitation, A = qcp-^ or q<pgjj, where cp-^ 
is the ionization potential and <p is the excitation potential (in volts). Since many 
levels may be excited, and the differential cross section should be different for each 

level. Therefore, the complete expression for the total momentum transfer cross sec- 
tion is given as follows: 
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where 



To determine a momentum transfer cross section, the differential cross section must be 
known. However, most experimental cross section measurements are of the total cross 
section, defined as 


cr,p(g) = jr^ sin x dx 


Experimental values of Orp , cr™, , and a™ as a function of electron energy for elec- 

^ef ^ex 

tron neutral collisions for several gases are presented in reference 7 (pp. 7, 102 and 111). 
Comparing the magnitudes of these quantities shows that, for electron energies ^ “g'^e) 
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lefis than 25 electron volts, a™, is at least an order of magnitude greater than Orp and 

CTrp . In reference 8 (p. 31), calculated values of cr are compared with Crr, . 

^ ex “^e£ e£ 

Since cr„ and cr™ are of the same order, it would appear that the same conclusion 

that holds for cr™ , or™ , and a™ vould also apply to , or^ , and , and as 
^e£ ^ex ”^ef "^ex 

a result a„ and ct will be neglected compared with cr_ in equation (B9). In 
”^I ”^ex ”^e£ 

this case then, equation (B4) becomes 



d'VggF F.ct (g) 

t s i m^j^ 


(BIO) 


The integration of equation (BIO) has been performed by several authors (ref. 23) for 
the case where Maxwellian velocity distributions are assumed for F^ and F^. In this 
case the result is 



cosh(^^^^ 

. <3^ 


e-gV^' dg 

(Bll) 

V/ 

2gw 





where w = Ug - u^. and /3^ = [(2kTg/mg) + (2kT^/m^)]. 

To carry out the g integration in equation (B 11), (t„ (g) must be known. Even if 

the integration is performed, the result will be a function of both w and j3, which is not 
a convenient form to be used in equation (B3); however, equation (Bll) can be simplified 
for the case w/jS « 1. For the Hall ion accelerator, this should be a good approxima- 
tion in the case of electron-neutral and electron- ion collisions: 
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For the case w/j3 — 0, equation (Bll) becomes 



2 

Now the integral will be a function only of the parameter ^ = (2kTg/mg) + (2kTj./mp . 
Equation (B12) can be written in terms of a momentum transfer collision frequency 
as follows ' 



(B13) 


where is defined in the following way: 
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(g)e 




dg 


(B14) 


Equation (B13) can be substituted into equation (B3) to obtain the equation of motion 
for the species s. The equation of motion for electrons in a plasma consisting of neu- 
trals, singly charged ions, and electrons is 
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m'g(Ug . V)Ug + + q(E + Ug X B) = - 


e o 


mm. 
e 1 


(iT - n) + 

^me+mo m ° + m. 


u . (H - u.) 
ei„ ' e r 


m 


(B15) 

where a steady state has been assumed, and V • ;^ = Vp^ since a Maxwellian velocity 
distribution has been assumed. 

For the case of electron neutral collisions (2kT /m ) + (2kT /m ) w (2kT /m ), and 

G G 0 0 GO 

thus the electron neutral momentum transfer collision frequency is 



(B16) 


In computing the integral in equation (B16), experimental values for (g) can be used. 


ef 


Since most of the experimental data is for cases where the neutral gas is initially at rest, 
the relative speed g is essentially the electron speed v . In these cases it is conven- 


ient to write equation (B15) in terms of the electron energy, (p =~ m v^ = — m g* in 


1 


volts: 


e e 



The total momentum transfer cross section for coulomb collisions can be computed 
analytically, provided that the lower limit in the expression 

= 27t /* (1 - cos x)cr. sin x dx 

“e£ Jo “ 


is cut off at some value, since the integral diverges at the lower limit in x* The value 
for when the integration is cut off at a lower limit corresponding to an impact pa- 

rameter equal to the Debye distance, defined in mks units as 
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1/2 



is given in reference 24 (p. 127) as 



= 4tt 


/ ^s^t 

\4ne^lJ.g^ 


2 


In A 


127re kT 

where e is the permittivity of free space and A = — rj^. Burgers assumes 

'Is'lt 

thermal equilibrium (T = T„ = T.) in obtaining this result; however, for the Hall ion ac- 

s r 

celerator the electrons have a much higher temperature than the ions and neutrals. To 
compute rj^ and In A, therefore, will be used in the expressions for r^ and A. 
Substituting the momentum transfer cross section in equation (B14) and assuming 
(2kTg/mg) » (2kT^/mj^) give the following results for the electron-ion momentum trans- 
fer collision frequency 


ei 


_ 4V^ / q^ 

' 3 \47re. 


m In A 




= 3.63x10 


-6 "e 

t3/2 


sec"^(mks) 


where equal charge densities have been assumed (n^ = n^) and 




1. 24x10"^ 



(mks) 


(B18) 
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APPENDIX C 


ELECTRON ENERGY EQUATION 


The energy equation for the species s can be obtained by multiplying the Boltzmann 

equation by ^ and integrating over v„ with the following result (ref. 24, eq. (5-43), 

z s s s 

p. 132) 


at 


m — 


3 Ps -s 


® ®\2 m 2 

s s 




+ 






'3 Ps .^s 


® ®\2 m 2 

s s 




+ V.,^g.Ug + V.Qs 


q n E • u = \ / ( — 

^sss s / I \ dr / 

/—j *'coii 


j (Cl) 


where Q, 




2— 3 

c c„F„ d v_ is the heat flow vector and p 
s s s s s 






s s 

is the scalar pressure. The term on the right side of equation (Cl) represents the energy 

gained by species s in collisions with all other species. 

Similar to the case for momentum transfer, the energy transfer term due to colli- 
sions can be written as follows: 


u 


<3F 
s 

^ 9t 


— mv^ d^v 


"coll 



Using equations (B5a), (B5b), and (B5c) gives the following results for (v^ 






= M V • (g - g*) - 


M 


(C3) 
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Using equation (C3) in equation (C2) and performing the x and e integrations results in , 




where o™ and a™ are the total ionization and excitation cross sections defined in 
^ i ^ ex 

appendix B. As was done in appendix B, the total momentum transfer cross section, 
will be approximated by the elastic total momentum transfer cross section. 

To perform the v and integrations a coordinate transformation will be made. 
Values for F and Tx must also be assumed, however. The most obvious approxima- 

SI ^ ^ 

tion for Fg and F^ is a Maxwellian about the mean velocities Ug and u^ at the species 
temperatures Tg and T^. Therefore, 


n 


^s = 


s 


^ 3 / 2,3 

7T ag 





V2a 


2 2 

where a = 2kT /m and a. = 2kTx/m.. The product of the distribution functions is 

S S S L L L 

transformed into the following expression 
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(g+w)^ 

2 .2 

S I o q 

TT^aga/ 


where the quantities appearing in this expression are defined as follows: 




C = V- 




w = Ug - Ut 




o2^2 
a 2 = Vt 

° o2 


M^2 M^2 

It is now possible to transform the integration in equation (C4) from Vg, to g 
since 


d^Vg d^v^ = 1 J| d^? d^g = dh d^g 

where the Jacobian, I j| = — = i. The Jacobian is calculated by using 

3(V,g) 0(C,D 

equations (B5) and the definition of f . Substituting the product of the distribution func- 
tions in equation (C4) and integrating over C result in 
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Changing the integrals to spherical coordinates and carrying out the angular integrations 
give 
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, Just as in the case for momentum transfer (appendix B), the limit w/)3 — 0 is taken. 
As a result the first integral vanishes: 




dg 


m^ 
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(l<Pl 


Jo Z— / Jo 


Q-is/13) 
ex 




(C7) 


Equation (C7) can be written in terms of the momentum transfer collision frequency, an 
ionization collision frequency, and an excitation collision frequency 


u- 


2 s s s 

•^coll 
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3m m. 

"si ^ 






(mg + m^y 


m + mg V 


(C8) 


where ^ is defined by equation (B14) and 
m 


. ^ r 

''"‘ex „V2^3 / 

^ex 




dg 0 -^, = 0 for g < g. 


(C9) 


2 2 

g^CTrj, (g)e~^ dg OP = 0 for g < 
ex ex 


(CIO) 


Substituting equation (C8) into equation (Cl) gives the energy equation for species s: 
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It should be noted that, since a Maxwellian velocity distribution function has been as- 
sumed, the heat flow vector vanishes and that V . • Ug)= ^Pg^g- Consider now the 

form equation (Cll) takes for electrons in a plasma made up of electrons, singly charged 
ions, and neutrals of the same species as the ions. For ionization and excitation of neu- 
trals by electron impact, equations (C9) and (CIO) in terms of the electron energy in volts, 
become 


“I ' V? Ue/ \2kTe 


/ oo 

^o^Tj' 


-q<p/kT 

{(p)(pe ® dcp (C12) 


and 




2/m ^V2 


"“ex 


^oO 


N Orp (<p)(pe 
ex 


-q<^/kT, 


d(p (C13) 


where; 


„ 1 ’^e_2 

<P = g 

2 q 


In deriving the preceding results, the condition, T /m » T /m has been used to 

o tJ fcJ L# U 

obtain j3 = 2kTg/mg. The ionization collision frequency for electron-neutral colli- 
sions is obtained by carrying out the indicated integration of Orp . A useful expression 

^I 

for is obtained at low temperatures by assuming Orp to be the following linear func- 
tion of (p ^ 

(No)std''Tj = - *^1^ 
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whfre aj is a constant, <pj is the ionization potential, and (No)g^^ is the number den- 

O 2 

sity at 0 C and 1 millimeter of mercury (133. 3 N/m ). Such a linear approximation for 

the ionization cross section is made in reference 6. Therefore, for an arbitrary number 

density (in reciprocal cubic meters) equation (C14) becomes 


Vt, ' - «’l> 


273 


(CIS) 


where is the neutral temperature in °K and p^ is the neutral pressure in newtons 
per square meter. In reference 11, von Engel gives values of aj for several gases when 
p is measured in millimeters of mercury and N or„ in reciprocal centimeters: hence, 

W O X jp 


aj = 


10' 


133. 3 N/m 


2 ^^I^v. Engel 


' ^ kg~^ ^ ®"Sel 


(C16) 


For argon. 


^^v. Engel = 


hence 


aj = 0 . 534 (A)(sec V(kg^)( m^) 

Substituting equation (Cl 5) in equation (C12) and integrating yield for small 

* q ' T„ ( y ( 2kT^y 


v3/2 . 

5.04xl0'^^ajNj^| e 


(kT. 




2kT, 


(C17) 


where 

Nq neutral number density, m 
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Tq neutral temperature, °K 

2 

Pq neutral pressure, N/m 

(pj ionization potential of gas (for argon, 15. 76 V), V 

To obtain a low temperature approximation for the electron neutral excitation colli- 
sion frequency, a linear variation of cr,p with electron energy will be made, just as 

^ex 

was done in the case of the ionization collision frequency. Only the first atomic level is 
considered. For argon, the data of Druyvesteyn and Penning (ref. 7, p. 102) were used 
to construct an expression similar to equation (C15). As a result, an expression of the 
same form as equation (C17) is obtained 


■'ex =5. 04X10-' 5 



-qcPexAT^ 


sec 


-1 


(C18) 


where 


^ex excitation potential of first level (for argon, 11.6V), V 
&ex defined the same way as aj in equation (Cl 4); for argon, 
^ex ^ (A)(secV(kg)(m^) 


Sovie and Dugan (ref. 25) have calculated the quantity q(^ji^j + X^<i^ex^ex)/^o 
argon, cesium, and helium as a function of electron temperature. To calculate v. 


for 

and 


^ex’ assumed a Maxwellian distribution for the electrons and used the method of 
Gryzinski (ref. 26) to calculate ionization and excitation cross sections. Their results 

are compared in figure 18 with values of 

^^‘^ex^ex ■*' ^I^lV^o argon calculated by 
using equations (C17) and (C18). Reasonable 
agreement is obtained for low values of the 
electron temperature. At an electron temper- 
ature of 14 electron volts, the results of refer- 
ence 25 and the calculated values begin to di- 
verge as a result of the assumed linear varia- 
tion of and 

0X 

in deriving 

14 electron volts, the calculated values are 
within 20 percent of the results in reference 25. 
For low temperatures, either evaluation of this 
inelastic loss parameter is equally justifiable, 
r- 10 , , .• .ion (X. iroxn oc At Mgh temperatures, the values calculated in 

Figure 18. - Inelastic power consumption for Argon as ® ^ ’ 

a function of electron temperature. reference 25 are more accurate. 



C7j with electron energy made 
Uj and However, below 
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